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Abstract — Large-scale power blackouts caused by cascading 
failure are inflicting enormous socioeconomic costs. We study the 
problem of cascading link failures in power networks modelled 
by random geometric graphs from a percolation-based viewpoint. 
To reflect the fact that links fail according to the amount of 
power flow going through them, we introduce a model where 
links fail according to a probability which depends on the number 
of neighboring links. We devise a mapping which maps links in 
a random geometric graph to nodes in a corresponding dual 
covering graph. This mapping enables us to obtain the first- 
known analytical conditions on the existence and non-existence 
of a large component of operational links after degree-dependent 
link failures. Next, we present a simple but descriptive model 
for cascading link failure, and use the degree-dependent link 
failure results to obtain the first-known analytical conditions on 
the existence and non-existence of cascading link failures. 

I. Introduction 

In August 2003, a massive power blackout rolled across 
northeastern U.S. and Canada, affecting some 55 million peo- 
ple. The blackout resulted from a three-hour-long cascading 
failure triggered by an undetected initial power line failure 
in Ohio. Later that year, 57 million Italians were left in the 
dark from a similar blackout. Around the world, the number 
of outages affecting large populations has steadily risen in 
the past decade, causing enormous economic losses JT|. As 
a consequence, the resilience and reliability of power systems 
have received increasing attention. Designing and constructing 
a smart, resilient, and self-healing power grid has become a 
crucial challenge for governments and industry. 

The essential dynamic of power blackouts is a process 
whereby the failure of a generating plant or a power line 
results in redistribution of the power flow load onto other 
nearby generators or power lines. If these other plants and lines 
then fail due to excessive load, then this process can further 
propagate and result in a cascading failure. A key question 
for the designer of the power grid is: can we predict whether 
a small number of initial failures will or will not trigger a 
cascading failure affecting the entire network? 

Previous work has examined the problem of cascading 
failure from a number of perspectives. In [2|, the authors use 
a linear model for power flows. They model transmission line 
failures as the result of overflow, and consider the dynamics of 
power networks in time. In [3|, the authors include the growth 
of power demand, engineering responses to system failures, 
and the upgrading of generator capacity. They characterize 
two types of blackouts, one which is due to transmission lines 



reaching load limits but incurring no line outages, the other 
which is due to multiple line outages. They use simulations 
to illustrate how their proposed model fit the data from 
American blackouts. In |4|. the authors characterize two kinds 
of critical operation points of the power system, due to the 
transmission line limit and the generator capacity limit, re- 
spectively. They find by simulation that operation near critical 
points could produce power law tails in the blackout size 
probability distribution. In [5 1, the authors propose a method to 
prevent cascading power failure by using integer programming 
techniques. They consider the cost of upgrading the capacity 
of transmission lines and calculate the optimal choice of 
link upgrade to prevent the expansion of failure. As seen 
by the brief review above, the current research on cascading 
failures in power networks is primarily based on simulation 
and optimization algorithms. Due to the complex nature of 
the dynamics involved in cascading failures, analytical results 
have been rather elusive. 

In 0, the authors propose a new approach to analyzing 
cascading failure based on the theory of percolation ll7l- lfTTI . 
The perspective adopted in [6| is as follows: in a power 
network subject to redistribution of load, the presence of 
cascading failure can be assessed by whether or not a small 
number of initial node (generator) or link (line) failures lead to 
a large connected component of failed nodes or links affecting 
the network globally. Percolation theory, which is concerned 
with assessing the global connectivity of networks, provides 
the appropriate viewpoint with which to make this assessment. 
The authors of |6| propose a simple but descriptive model for 
cascading node failure, based on node susceptibility thresholds 
and degree-dependent node interactions. Within the context 
of random geometric graphs, |6| presents the first analytical 
conditions for the existence and non-existence of cascading 
node failures. 

In this paper, we adopt the viewpoint presented in (6). 
Instead of examining cascading failure of nodes, however, 
we focus on cascading link failures. This is motivated by the 
fact that blackouts in power grids are often triggered by link 
failures (line faults), which cause power flow redistribution 
onto neighboring links, which when overburdened, can fail 
and propagate the cascade. As in (6), we shall focus on 
networks modelled by random geometric graphs, where nodes 
are distributed spatially according to a Poisson point process 
with constant density A, and any two nodes within given 
distance r of each other share a link. In doing so, we are 



motivated by the following consideration. First, the topology 
of actual power grids is not readily available. Second, we use a 
random graph to model the ensemble topology of power grids. 
Third, there is evidence that geometry plays an important role 
in the topology of electrical power grids lfT2l . especially with 
respect to cascading dynamics. Finally, our work is intended 
to provide insight into the dynamics of cascades within large- 
scale power networks, and random geometric graphs provide 
a concrete setting in which analytical results may be obtained. 

This paper proceeds as follows. After presenting some pre- 
liminary results on random geometric graphs and continuum 
percolation, we study random geometric graphs in which 
each given link fails (independently) with a probability which 
depends on the number of neighboring links (links which share 
an end node with the given link). This models the phenomenon 
in power networks where the probability of a link failure 
increases with its load in terms of power flow. The larger 
the number of other links which share an end node with a 
given link, the more power flow is likely to traverse the given 
link, and therefore the higher the probability that the given 
link fails. We present the first-known analytical conditions 
on the existence and non-existence of a large component 
of operational links after degree-dependent link failures. In 
order to establish these results, we devise a new map which 
maps links in a random geometric graph to corresponding 
nodes in a dual covering graph. Finally, we present a simple 
but descriptive model for cascading link failure, and use the 
degree-dependent link failure results to obtain the first-known 
analytical conditions on the existence and non-existence of 
cascading link failures. 

II. Network Model 

We use a random geometric graph to model the power 
network. A rigorous model is as follows. Let H\ be a 
homogeneous Poisson point process with density A > in M 2 . 
Let G(H\,r) be the (infinite) undirected graph with vertex set 
"Hx and with undirected links connecting all pairs {Xj,Xj} 
with ||Xj - Xj-|| < r (r > 0) where X;,Xj e U\. Due to 
the scaling property of random geometric graphs |9), in the 
following, we focus on G(Hx, 1)- 

A fundamental result for random geometric graphs concerns 
a phase transition effect whereby the macroscopic behavior of 
the network is very different for densities below and above 
some critical value A c . For A < A c (subcritical or non- 
percolated), the connected component containing the origin 
contains a finite number of points almost surely. For A > A c 
(supercritical or percolated), the connected component con- 
taining the origin contains an infinite number of points with a 
positive probability [7|-[9]. 

Let Wx.o = Hx U {0}, i.e., the union of the origin and 
the homogeneous Poisson point process with density A0 As 
discussed, a phase transition takes place at the critical density. 

1 Note that in a random geometric graph induced by a homogeneous Poisson 
point process, the choice of the origin can be arbitrary. 



Definition 1. For G(Hx,o,^), the percolation probability 

Poo (A) is the probability that the connected component con- 
taining the origin has an infinite number of nodes of the graph. 
The critical density A c is defined as 

A c = mf{A>0: Poo (A) > 0}. (1) 

It is known that if A > A c , then there exists a unique infinite 
component in G(Hx, 1)- A fundamental result of continuum 
percolation states that < A c < oo 1T31 . Exact values of 
A c and Poo(A) are not yet known. Simulation studies show 
that 1.43 < A c < 1.44 

We have defined the random geometric graph in the infinite 
graph setting. The finite analog of G(Hx, 1) can be defined 
as follows |9). Let X n = {Xi,...,X n } be n points inde- 
pendently and uniformly distributed in a square of area n/X 
in M 2 . Form a graph G(X n , 1) by placing a link between 
any two nodes i and j such that ||Xj — Xj|| < 1. Let 
L±(G(X n , 1)) be the size of the largest connected component 
in G(X n , 1). Then it can be shown that as n —> oo with A 
fixed, n~ 1 Li(G(X n , 1)) — > Poo(A) in probability. Thus, if 
A > A c , there exists a unique largest component of size 0(n) 
in G(X n , 1). On the other hand, it can be shown that if A < A c , 
then the largest component in G(X n , 1) can only have size 
O(logra) (9). 

While finite random geometric graphs are clearly of more 
interest for practical applications, it is more convenient to 
consider infinite graphs for expositional purposes. For this 
reason, we will concentrate on infinite graphs in the rest of 
the paper, while keeping in mind the correspondence between 
the finite and infinite graphs. 

III. Degree-dependent Random Link Failures 

In the power grid, large-scale blackouts are often triggered 
by power line faults. This was the case in the massive 2003 
North American blackout. In our context, line faults can be 
modeled as link failures. Thus, in order to further understand 
blackouts in the power grid, it is important to characterize the 
resilience of large-scale networks to link failures. 

In power networks, the probability of link failure typically 
increases with its load in terms of power flow. The larger 
the number of other links which share an end node with a 
given link, the more power flow is likely to traverse the given 
link, and therefore the higher the probability that the given 
link fails. To model this effect, we consider degree-dependent 
random link failures. We first define the degree of a link. 

Definition 2. The degree of link denoted by d(i,j), is 

the number of other links which share an end vertex with 
That is, 

d(i, j) = di + dj — 2, 
where di and dj are the node degrees of i and j. 

Consider a scenario where each link in G(HxA) fails 
(independently) with a probability q(k) that depends on the 
link's degree k. Let G("Ha, 1, <?(■)) be the remaining graph 



after degree-dependent link failures. In order to study the con- 
nectivity properties of G(1~L\, 1, <?(•)), we construct a mapping 
where the bond percolation process in G(H\, 1) is mapped to 
a site percolation process in a covering graph G C (W\, 1). 
The covering graph of a graph G is defined as follows. 

Definition 3. G C (V,E) is called the covering graph of 
G(V,E)if 

1) Each node in G c denotes a link in G, and each link in 
G is denoted by a unique node in G c - 

2) Two nodes in G c share a link if and only if the cor- 
responding links in G share a common end vertex (in 
G). 

Comparing Definitions [2] and [3] we have the convenient 
fact that the degree of a link in G is equal to the degree of its 
corresponding node in G c . 

It should be clear that the covering graph of a given graph is 
not unique. Moreover, in the context of geometric graphs, most 
covering graphs are of little importance, for they hold only the 
connection information, losing all the geometric information. 
In order to utilize the geometric information in the original 
graph, we need to place additional constraints on the geometric 
location of nodes in the covering graph. Fortunately, we find a 
construction for the covering graph which is particularly useful 
for random geometric graphs with random link failure. In this 
construction, the location of nodes in the covering graph is 
directly determined by the node locations in the original graph, 
thereby maintaining the geometric information. 

Given G(H\,l), we construct the covering graph 
G C (H\, 1) as follows. 

1) Consider the circles with radius 1/2 and the nodes in 
G(H\, 1) as the centers. 

2) Let the mid-points of each link in G (H\ , 1) be the nodes 
in G c {Hx,l). 

3) Two nodes in G C (H\, 1) share a link if and only if they 
lie in the same circle of radius 1/2 in G(H\, 1). 

The construction is illustrated in Figure Q] 




Fig. 1. Construction of the covering graph: the red nodes and links belong 
to G(H\, 1); the blue nodes and links belong to G C (H\, 1)- 

Theorem 1. The graph G C {%\,1) constructed above is a 
covering graph of G(H\, 1). 

Proof: First, each link in G(H\,1) corresponds to a 
unique node in G C (H\,1), which is the link's midpoint. 



Second, since the length of links in G(H\, 1) is less than 
or equal to 1, two nodes in G C {H\, 1) lie in the same circle 
(of radius 1/2) in G(H\, 1) if and only if their corresponding 
links in G(W\,1) share an end vertex. Thus G C (H\,1) is a 
covering graph of G(Ji \, 1). ■ 
One important consequence of the above construction is that 
the nodes of the covering graph G C (H\, 1) is also driven by 
a Poisson point process, albeit of a different density. 

Theorem 2. The nodes ofG c {'H\ 1 1) are distributed according 
to a Poisson point process with density X' = 

Proof: To prove that the nodes in G C (H\, 1) are driven 
by a Poisson process, it is most convenient to focus on 
the finite random geometric graph G(X n ,l) where X n = 
{Xi,...,X„} are n points independently and uniformly 
distributed in a square A of area n/X in R 2 . Now consider 
the covering graph G c (X n , 1) of G(X n ,l) obtained by the 
construction above. For any i = 1, . . . , n, conditioned on the 
location of Xj, the midpoint of any link between X^ and Xj 
in G(X n , 1) is uniformly distributed in the circle centered at 
Xi with radius 1/2. Since the X;'s are independently and 
uniformly distributed in A, so are the midpoints, which are 
the same as the nodes of the covering graph G c (X n , 1). 

To calculate the density of G c (X n , 1), we note that the mean 
degree of any node X^ in G(X n , 1) (ignoring border effects 
which diminish as n tends to infinity) is irX. The summation 
of the mean degrees of all nodes Xi,X2, ...,X n equals nirX. 
Now the number of edges of the graph should equal 1/2 of 
the sumation of the degrees of all the nodes. Thus the mean 
number of edges in G(X n ,l) equals nirX/2. Therefore, the 
density of the links in G(X n , 1), which is also the density of 
the nodes in the covering graph G c (X n , 1), is given by 

, nirX/2 7rA 2 

The theorem now follows by taking the limit as n —> oo with 
A fixed. ■ 
It is important to remark that even though the nodes of 
the covering graph G C (H\, 1) are distributed according to a 
Poisson process, G C (H\, 1) is not a random geometric graph. 
To see this, note that two nodes in G C (H\,1) being within 
distance 1 of each other does not imply that there exists a link 
between the nodes. Nevertheless, G C (H\, 1) is a subgraph of 
a random geometric graph with the same node locations as 
G c (Hx,l). 

We now derive analytical conditions on the existence and 
non-existence of an infinite connected component of oper- 
ational nodes in G(H\,l,q(-)) after degree-dependent link 
failures in of G(H\,1). The key techniques we use are the 
mapping from a bond percolation process in the original graph 
to a site percolation in the covering graph, and that link failures 
affect network connectivity less than node failures at the same 
level. 

Our first main result gives a sufficient condition for there to 
exist an infinite connected component of operational nodes 



in G(H\, 1, (?(•)) after each link in G("Ha,1) fails (inde- 
pendently) according to a degree-dependent probability q(k), 
where k is the degree of the link, as given in Definition [2] 

Theorem 3. For any Ai > A c and G(H\, 1) with A > Ai, 
there exists k\ < oo which depends on A and Ai, smc// f/iflf if 

?(*) < 1 - ^r, far alll<k< ki, (2) 
A 

then with probability 1, there exists an infinite connected 
component in G(H\, 1, #(•)). 

An interesting implication of Theorem [3] is that even if all 
links with degree larger than k\ fail with probability 1, an 
infinite component still exists in the remaining graph as long 
as (fJJ is satisfied. 

Proof: The proof uses the concept of the covering graph, 
the relationship between bond percolation and site percolation, 
as well as key ideas from the proof of Theorem l(i) in J6). 
Due to space constraints, we shall only outline the key steps. 

Consider two random failure models in G(W\,1). In the 
first model, each node fails (with all associated links) indepen- 
dently with probability 1-^-. Let Gf te {Hx, 1) be the remain- 
ing graph. In the second model, each link fails independently 
with probability 1 — 4^. Let Gi(H\, 1) be the remaining graph. 
Since Ai > A c , G± te {Jl\, 1) is a random geometric graph in 
the supercritical regime. Let G c (%\, 1, q(-)) and Gi c (%\, 1) 
be the covering graphs (according to the construction given 
above) of G(H\, 1, <?(•)) an d G\(H\, 1), respectively. 

In the following, we use the square lattice construction 
from [6|. Please refer to |6| for further details on the con- 
struction. Consider the square lattice £ = d ■ Z 2 , where d 
is the edge length. As in [6|, define event A s ^ te (d) for each 
horizontal edge a in £ as the set of outcomes for which the 
rectangle R a is crossed from left to right by a connected 
component in Gf te (H\,l). If A^ lte (d) occurs, we say that 
rectangle R a is a good rectangle, and edge a is a good edge. 
Let pf e (d) = Pr {Af e (d)). Define Af e {d) similarly for all 
vertical edges by rotating the rectangle by 90°. 

Further define event B^ ite (d) for edge a in £ as the 
set of outcomes for which (1) A^ te (d) occurs, and (2) the 
left square S~ and the right square 5+ are both are both 
crossed from top to bottom by connected components in 
Gf lte (W\,l). If B^ te (d) occurs, we say that rectangle R a 
is a complete rectangle, and edge a is a complete edge. Let 
pf e (d) = Pr(B s a lte (d)). Define B s a lte (d) similarly for all 
vertical edges by rotating the rectangle by 90°. An example 
of a complete rectangle is illustrated in Figure [2] 

We further define complete events {B' a (d)} with respect to 
Gi(H\,l) in the same way as we defined complete events 
{B s a tte (d)} with respect to Gf e (H\,l). Similarly, define 

2 Here, a rectangle R = [0:1,2:2] x [2/1,2/2] being crossed from left to 
right by a connected component in Gf Ue ('H\,l) means that there exists 
a sequence of nodes Vi,V2, ■■■,v m £ G\ lte CH\,1) contained in R, with 
|x„ ; - x„ i+1 1 < l,i = 1, ...,m - 1, and < x(«i) - xi < 1/2,0 < 
X2 — x(v m ) < 1/2, where x(tli) and x(v m ) are the ^-coordinates of nodes 
vi and v m , respectively. A rectangle being crossed from top to bottom is 
defined analogously. 



Fig. 2. An example of a complete rectangle. The circles have radius 1/2. 

complete events {B"(d)} with respect to Gi c (%a, 1). the 
covering graph of Gi^a,!), and define complete events 
{B a (d)} with respect to G C {U\, 1, <?(•))• 

We shall now carry out a series of stochastic couplings to 
obtain our result. Since the critical link failure probability (for 
i.i.d. link failures) is greater than or equal to critical node 
failure probability (for i.i.d. node failures) in G("Ha,1), we 
can stochastically couple Gi(H\,l) and Gf te (H\,l) such 
that the existence of crossings defined in events {B^ te (d)} 
for G s i te (H\, 1) implies the existence of crossings defined in 
events {B' a (d)} for G\(T-L\, 1). Next, by using the construction 
of the covering graph Gi c (W\, 1), it is straightforward to ver- 
ify that the existence of crossings defined in events {B' a (d)} 
for Gi('Ha,1) implies the existence of crossings defined in 
events {B'^d)} for G lc (Hx, 1)0 

Let p"(d) = Pr(B"(d)). By the coupling argument above, 
we have p"(d) > p s ^ te (d). From (6), we know that pf te (d) 
converges to 1 as d —> 00 when Gf lte ('HA,l) is in the 
supercritical phase. Thus, p"(d) also converges to 1 as d — > 00. 

Define 

d(Y, Ai) 4 inf Id > 4 : p'^d) ~ (d +2 ) (*» + 2 ) A' > 1 ~ 9 °\ ' 

(3) 

where A' is the density of G C (T-L\, 1), the covering graph of 
G(W\,1), and qo = g+ ^/^ - Now choose the edge length of 
C as d = d(X', Ai). 

Define event C a (d) for each horizontal edge a in £ as the 
set of outcomes for which the number of nodes of G C (H\, 1) 
in R' a is strictly less than 

fcl 4 2 (:aw +2 )(M^) +2 )y, (4) 

where R' a is the rectangle R a extended by 1 in all directions. 
Define C a (d) similarly for all vertical edges by rotating the 
rectangle by 90°. If C a (d) occurs, we call rectangle R a and 
edge a efficient. Let p e (d) = Pr(G a (d)). 

We say an edge a in £ is open if and only if it is both 
complete with respect to Gi c ("Ha, 1) and efficient with respect 
to G c ("Ha, 1), i-e. when events B"(d) and C a {d) both occur. 
An edge is closed otherwise. 

3 To be absolutely precise, the definition of a crossing needs to be altered 
slightly here, requiring the x'-coordinates of the first and last node in the 
rectangle to be within 1 of the boundary, rather than 1/2. However, this 
alteration does not affect our results below. 



When C a (d) occurs for edge a in C, no node of 
G c (%\, 1, (?(•)) in R a has degree strictly greater than k\ in 
G C (H\, 1). In addition, if g(fc) satisfies (fJJ, a link in G(H\, 1) 
with degree fc, 1 < k < k% (or equivalently a node in 
G c ('Ha,1) with degree k, 1 < fc < fci), survives with a 
probability greater than or equal to ^ in the degree-dependent 
failures model. On the other hand, for the independent random 
failures model, a link in G(W\,1) (or equivalently a node 
in GciTix,!)) survives with probability exactly equal to 4^. 
Thus we can couple G C (H\, 1, q(-)) with Gi c (W\, 1) so that 
the existence of crossings defined in events {B"(d)} for 
G\ C {T-L\, 1) implies the existence of crossings defined in events 
{B a (d)} for G C {H\, 1, <?(■))■ Therefore, a path of open edges 
in C implies a connected component in G C (H\, 1, <?(•)). 
As in [6], we find 



and 



Po {d) = Pr(K(d) n c a ) > p'^d) +p e {d) - 1, 



Pe(d) > 1 



(5) 



(6) 



(f^M +2 ) (M^il + 2 ) A' 
By (0, © and ©, we have 

Po(d) > p c (d)+ Pe (d) -l>l-q . (7) 

Finally, we use the dual lattice technique in |6| to prove 
the existence of an infinite open edge cluster in C, which 
implies the existence of an infinite connected component in 
G c i1-L\, 1, ?(•))> an d therefore an infinite connected compo- 
nent in G(H\, 1, £?(•))• ■ 

Theorem[3]provides a sufficient condition for G(H\, 1, q(-)) 
to have an infinite component. The next theorem provides 
a sufficient condition for G(H\,l,q(-)) to have no infi- 
nite component. Thus, it provides a necessary condition for 
G(W\, 1, q(-)) to have an infinite component. Here again, the 
concept of the covering graph allows us to leverage Theorem 
l(ii) of |o| to prove our new result. 

Theorem 4. Given G{U X , 1) with X > X c , let X' = irX 2 /2 be 
the density of the covering graph G C {T-L\, 1). If either 
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when q(k) is non-decreasing in k, or if 
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(l-q(m + k-l) k )j < - (9) 

when q(k) is non-increasing in k, then with probability 1, there 
is no infinite connected component in G(H\, 1, (?(•)). 

Proof: The existence of an infinite connected compo- 
nent in G(H\, 1, <?(•)) a f ter degree-dependent link failures in 
G(W\,1) according to failure probability q(k) (where k is 
the link degree) is equivalent to the existence of an infinite 



connected component in the covering graph G C (H\, l,q(k)) 
after degree-dependent node failures in G C (H\, 1) according 
to failiure probability q(k) (where k is the node degree). 

In G c (T-l\, 1), the nodes are Poisson distributed with density 
A'. The length of any link is at most 1, and any two nodes with 
distance strictly larger than 1 do not share a link (although 
two nodes with distance less than 1 may not necessarily 
share a link). Now consider the random geometric graph 
G r (H\>, 1) with density A' where the node locations are the 
same as those for G C (H\,1). Let G r (H\> ,l,q(-)) denote 
the remaining graph after degree-dependent node failures 
in G r (H\>,l) according to failure probability q(k). Since 
G C (H\,1) is a subgraph of G r (H\> ,1), we can stochasti- 
cally couple G r (H\,l,q(-)) and G C (H\, 1, «?(•)) such that 
if there does not exists an infinite connected component in 
G r (W\, 1, q(-)), then there also does not exist an infinite 
connected component in G C (H\, 1, </(•)), By Theorem l(ii) of 
[6|, if either ^ or © is satisfied, then G r (H\, 1, <?(•)) does 
not have an infinite component, and therefore G C (H\, 1, <?(•)) 
does not have an infinite component either. ■ 

iv. Cascading Link Failure 

Cascading failure (blackouts) in power networks often result 
from a small number of initial line faults triggering many 
more line failures affecting the entire network. To understand 
the underlying dynamics of such phenomena, we focus on 
cascading link failures in large-scale networks modelled by 
random geometric graphs. 

Consider a network modelled by a random geometric graph 
G(W\,1) with A > A c , where an initial failure seed is 
represented by a single failed link. We are interested in 
whether this initial small shock can lead to a global cascade 
of link failures, defined as follows. 

Definition 4. A cascading link failure is an ordered sequence 
of link failures triggered by an initial failure seed resulting in 
an infinite component of failed links in the networks. 

To describe cascading failures, we use the following simple 
but descriptive model similar to the model proposed in for 
cascading node failures. Assume that there is a "susceptibility 
threshold" ipi 6 [0, 1] associated with each link I, where 
the ipi's are i.i.d. random variables with probability density 
function f(ip). Due to local redistribution of power flow load 
after link failures, each link / fails if a given fraction ipi of its 
neighboring links (the links which sharing an end node with 
T) have failed. The order of the failure sequence is then the 
topological order determined by the location of the initial link 
failure and the threshold i/ji of each node /. 

Note that the initial link failure can grow only when some 
neighboring link, say m, of the initial failure seed has a 
threshold satisfying ip m < where k m > 1 is the link 
degree of m. Such a link is called vulnerable. The probability 
of a link being vulnerable is 



(10) 



where F^(-) is the cumulative distribution function of ipi. 
When the initial failure seed is directly connected to a com- 
ponent of vulnerable links, all links in this component fail. 
Thus, a cascade of failed links forms when the network has 
an infinite component of vulnerable links and the initial failure 
seed is either inside this component or adjacent to (share a link 
with) some link in this component. 

On the other hand, if link / has a threshold satisfying ipi > 
where k\ is the link degree of I, then link I will not fail 
as long as at least one neighboring link is operational. Such a 
link is called reliable. Otherwise, if ipi < j^ 1 , we call link I 
unreliable. For k > 1, the probability of a node being reliable 
is given by 

a k = 1 - F 4 , (^—^-) = I (ID 



Note that when two reliable links share an end node and 
neither is an initial failure seed, no matter what else happens 
in the network, they remain operational. 

The following two theorems present our main results on 
cascading link failure in random geometric graphs. The results 
depend crucially on the degree-dependent link failure results 
in Theorems [3] and |U 

Theorem 5. For any \\ > X c and G(W\, 1) with A > Ai, 
there exists k\ < oo depending on A and Ai such that if 



(12) 



then with probability 1, there exists an infinite component 
of vulnerable links in G(H\,1). Moreover, if the initial link 
failure is inside this component or adjacent to some link in 
this component, then with probability 1, there is a cascading 
link failure in G(H\, 1). 

Proof: We view the problem as a degree-dependent link 
failure problem where a vulnerable link is considered "oper- 
ational" and a non-vulnerable link is considered a "failure." 
In this model, each link with degree k fails with probability 
1 — pk- Applying Theorem [3] directly, we obtain the result. ■ 

Theorem 6. For any G{T-L\, 1) with A > A c , if 
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where A' = 7rA 2 /2 and F^,(—oo) = by convention, then with 
probability 1, there is no infinite component of unreliable links. 
As a consequence, with probability 1, there is no cascading 
link failure in G(H\, 1) no matter where the initial link failure 
is. 

Proof: We apply Theorem [4] Regard an unreliable link 
as "operational" and a reliable link as a "failure". Then, 
<jfc becomes the failure probability q(k) in the context of 



Theorem |4] Since ak is non-increasing in k, we replace 
q{m + k - 1) in © with a m+k ^ = 1 - F^ (^Ef) 
and obtain dT3l . Thus, by Theorem @1 when ( TT3l holds, with 
probability 1, there is no infinite component of unreliable links 
in the network. Next, using the covering graph and techniques 
from the proof of Theorem 2(ii) in [6 |, we show that if there is 
no infinite component of unreliable links, then with probability 
1 , there is no cascading link failure no matter where the initial 
link failure is. ■ 

v. Conclusion 

In this paper, we studied the problem of cascading link 
failures in the power grid from a percolation-based viewpoint. 
To reflect the fact that links fail according to the amount 
of power flow going through them, we introduced a model 
where links fail according to a probability which depends on 
the number of neighboring links. We introduced a mapping 
which maps links in a random geometric graph to nodes in a 
corresponding dual covering graph. This mapping enabled us 
to obtain the first-known analytical conditions on the existence 
and non-existence of a large component of operational links 
after degree-dependent link failures. Next, we presented a 
simple but descriptive model for cascading link failure, and 
used the degree-dependent link failure results to obtain the 
first-known analytical conditions on the existence and non- 
existence of cascading link failures. In particular, we revealed 
the important roles played by vulnerable and reliable links. 
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